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Abstract 

We investigate conditions under which multiatom absorption of a 
single photon leads to cooperative decay. Our analysis reveals the 
symmetry properties of the multiatom Dicke states underlying the 
cooperative decay dynamics and their spatio-temporal manifestations, 
particularly, the forward-directed spontaneous emission investigated 
by Scully et al. 

Dicke pioneered the notion of cooperative spontaneous emission by a col- 
lection of atoms, highlighted by the "superradiant" A^^-scaling of the 
emission rate into resonant modes fl]. His work has prompted numerous 
studies of the dependence of cooperative spontaneous emission upon the 
initial preparation (cooperative dipole moment and excitation), as well as 
the spatial symmetry and density (interatomic distances) of the multiatom 
sample [21 [H HJ [5l El [7] . These factors determine the degree of cooperativity, 
which may be attributed to multiatom interference of radiated photons. This 
cooperativity may range from maximal enhancement (superradiance or su- 
perfluorescence), first observed in p], to maximal suppression (subradiance) 
[H El [7], corresponding to constructive or destructive interference, respec- 
tively. Similar effects in neutron scattering on crystalline lattices have also 
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been predicted [S]. The most powerful systematic means of classifying such 
effects is the multiatom Dicke-states basis, which embodies the symmetry 
properties of the system [6l [7]. The difficulty is the large multiplicity of 
Dicke states for iV 3> 1 and their contrasting symmetry properties. The 
cooperative characteristics of the emission may be obscured if many Dicke 
states of different symmetry become entangled or mixed by the dynamics of 
the process. Hence, the recent analysis of Scully et al. [9], implying that 
atoms sharing one photon absorbed at t = should subsequently re-emit 
it only in the direction of the absorbed photon (the forward direction) is a 
nontrivial and perhaps counterintuitive manifestation of cooperativity in a 
large volume of randomly distributed atoms without an initial cooperative 
dipole moment. This aspect is implicit in the far- field analysis of Rehler and 
Eberly [2l[T0]. 

Here we wish to revisit the foregoing problem using the methods of group 
theory [HI [12], so as to elucidate the following important questions: Which 
Dicke states become populated or excited during the process? What are 
their symmetry properties and how are they reflected in the spatio-temporal 
buildup of cooperative spontaneous emission? 

Consider the absorption of a single photon by a collection of N atoms, 
initially in their internal ground state. The photon frequency cko is chosen to 
be exactly at resonance with the atomic transition frequency ujq. We assume 
that the spatial distribution of the atoms is spherically symmetric. For the 
sake of obtaining analytic results, we use a Gaussian probability distribution 
of the random atomic positions Yj inside the sphere: 

V{v) = {^R,)-^eM-rVRl)- (1) 

Here Rq stands for the typical size of the spherical sample. 

We expand the wave function of the system "field + matter" as follows: 

N 

\m) = Y: a,e-^-"%) + Y: «:ke-'^'*|k), (2) 
i=i k 

where \ej) denotes all the atoms in the ground state, except for the excited jth 
atom, the electromagnetic field being in the vacuum state, and |k) denotes 
all the atoms in the ground state, with one photon present in the k mode. 
In the present treatment we neglect the effects of photon polarization. The 
Schrodinger equation then reads 

iaj = ^^kC^'^^'^Kk, (3) 

k 
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N 

ik^ = {ck-uJo)Kk + J29le-''^'''aj, (4) 

i=i 

where (7k is the atom-photon couphng constant for the k mode. The initial 
amphtudes in ([2]) are expressed by 

a,(0) = M0) = 0. (6) 



The integration of (jl]) permits us to express through a time integral in- 
volving aj. Then, pulling aj out of the time integral on the assumption that 
aj decays slowly on the scale of the cooperation time (see below) , we convert 
([3]) into the following equation: 



N 



where 



sin(co'o — ck)t .cos(ct;o ~ ck)t — 1 
ujq — ck ooq — ck 



(7) 



The second term in the square brackets in ([7]) diverges and thus requires 
renormalization. It corresponds to the Lamb shift of the optical transition, 
which is, in principle, different for each of the collective excited states 
given by the mutually orthogonal linear combinations of \ej). This varying 
part of the Lamb shift is of co-operative origin O [6] . It can be calculated 
upon noting that the finite sample size and the corresponding wave vector 
spread Ak ~ Rq^ remove the ultraviolet divergence characteristic of the 
(unrenormalized) Lamb shift for a single atom [TH]. We then find that the 
contribution of co-operative effects to the Lamb shift of a collective state is 
approximately fco-Ro times smaller than the decay rate of this state. Hence, 
the variance of the Lamb shift is insignificant compared to the collective 
decay rate. It allows us in what follows to include, as usual, the Lamb shift 
into the definition of the transition frequency c^o, retaining only the real 
''™i^-"cfc^^* term in the square brackets on the right-hand-side of (I7j). 

For t ^ uJq^ we can then substitute in ([7]) ~ it6{ujo — ck) and 

the expansion k ^ n{kQ + 6k), where Sk -C ko and n is the unit vector in the 
k direction, thereby obtaining the following estimation: 

^n'{t) ^ J ^vr|(7krf^(cA;)|k=„fcoe^'""(^^-'^')e[ct - |n(r, - r,0|], (8) 
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where Q{ck) is the density of photon modes, 9(x) is the Heavyside step 
function, which is equal to 1 for x > 0, | for x = 0, and for x < 0. 

One can see from ([8]) that if ct is much smaller than the mean interparticle 
distance then the decay matrix Tjj/ is diagonal, all its diagonal elements being 
the same, equal to the single-atom decay rate 



y ^7r|^kr^(cfc)|k=nfco, (9) 
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indicating the total absence of cooperativity. As the time t increases, co- 
operativity is established among increasingly more atoms. At t ~ i?o/c the 
collective regime of radiation is fully established. In what follows we consider 
Rq/c. In this limit Eqs. ([6], [7]) reduce to 



N 



a 



j'=i 



sin ko 


- r 











71 E 77'^ r cyr, t » Ro/c. (10) 



It is extremely difficult to calculate the exact eigenvalues and eigenstates 
of ffTOj) . Therefore, in what follows we introduce states that closely approxi- 
mate the eigenstates of the problem. To this end, it is convenient to include 
the phase factors associated with the incident photon momentum into the 
definition of the excited states, 

e^^ne,) - |e,). (11) 

The corresponding new probability amplitudes are 

^ e-^^o^a,. (12) 

In this basis, (JTOll takes the form 

TV 

Pj = -71 E Fir, - r,0/3/, (13) 
i'=i 

where 

ko\rj - rj>\ 

The powerful theory of the permutation group representations and their 
characterization by Young tableaus [H] provides a general recipe for con- 
structing the functions of arbitrary symmetry with respect to permutations. 
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Our case is relatively simple, since we deal with a system of two-level atoms. 
The relevant Young tableaus contain one of two rows only, i.e., they are de- 
noted, respectively, by {A^} or {A^ - A^', A^'} with A^ - A^' > A^' > 0. To 
construct corresponding wavef unctions, we apply the method described in 
[12] . Consider an operator W defined as follows. If VV is apphed to the wave 
function of A^ atoms in their internal ground state, the result is YljLi kj). If 
the jith, j2th, ... , jmth atoms are initially excited, then 

W\ej,)\ej^) . . . \ejj = \ej,)\ej^) . . . |ej„)^'|ei), (15) 

j 

where J^/ denotes the sum over j ^ ji, j2, . . . , jm- This VV commutes 
with any product of generalized permutation operators, which interchange 
not only the internal-state atomic variables, but the coordinate-dependent 
phase factors as well: 

f le;), i = j 
Ojilci) = I \ej), 1 = 1 . (16) 
I |ei), i^ 

Thus VV conserves the symmetry type of the state. Upon applying VV to the 
wave function of A^ atoms in the ground state, which is totally symmetric 
with respect to generalized permutations, we obtain the totally symmetric 
state with one atom excited 

1 ^ 

V A^ j=i 

which constitutes a one-dimensional group representation characterized by 
the Young tableau {A^}. We now construct from A^ linearly independent 
states |ej), by orthogonalization to the A^ — 1 states 

e,) 

(18) 

which comprise the basis of the irreducible representation characterized by 
the Young tableau {A^ — 1, 1}. The states flTSl) are normalized to 1 and 
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N-l 



) 



1 + (l/VAT) 
N-l 



= J2fh)^ 1 = 1,2,..., N-l, 



orthogonal to each other and to 10''^^^). Any product of pairwise operators 
f ll6p transforms any of the wave functions (1181) into a hnear combination 
of these functions, without adding terms containing Icp^^^). This procedure 
may be extended to the construction of doubly-excited states with the Young 
tableau {A^ — 2, 2} and so on. 

Expanding the wave function of the singly-excited atomic states as 

N-l 

|^e..)=cW|0^^>)+E^r'''l</>P"''^ (19) 

1=1 

we arrive at the following set of equations, whose terms are explained below: 

7V-1 



C 



= -7co.c^^> - E ^r4'^"''^ (20) 



1=1 

N-l 



c. 



{N-1,1} ^ _„JN} _ n,„r^^-i'i> 



..ci^> - E Qu'c}'-'^'^, (21) 



l'=l 

satisfying the initial conditions 

c^^>(0) = l, cp"^'^>(0) = 0, / = 1, 2, ... AT-i. (22) 
The collective decay rate of the fully symmetric state is then found to be [cf. 

m] 

^ N ^ 

7co. = iEE^(r.-r,'). (23) 

The coupling term in (1201) . mixing the fully symmetric state and the Ith state 
of lower symmetry, appears because the states f|T7|) . f|T8l) are not the exact 
eigenstates of the decay operator on the right-hand-side of f|T3|) . However, 
they provide a good approximation thereof, because the mixing is character- 
ized by the coupling strength 

^^ = 7=EE/j^(r.-r,0, (24) 

and is weak (see below) compared to the collective decay rate. The coeffi- 
cients 

N N 

Qiv = liT.i: f^ffFi^, - r.') (25) 
j=ij'=i 



{N} {N-l,l} 
• ••T1 gTT- --T1 




ground> 



Figure 1: Schematic representation of the processes described by (fTTl [TSl) . 
Young tableaus characterizing symmetry of the satates are indicated. 

in fl2T|) describe the decay of the Ith state of the {A^ — 1, 1} symmetry type 
if I = I', or the mixing of states with I ^ I'. The level scheme, the coupling 
and decay channels described by fl2Ul I^T]) are shown in Figure 1. 

Let us first calculate the collective decay rate 'jcoi- Since N ^ 1, perform- 
ing the double sum in (l23l) and dividing it by is equivalent to averaging 
over the atomic positions r^, rjr. 

Icoi = liN J dhj J dhy V{v^)V{vy)F{v^ - V,,). (26) 

Here we assume that there is no correlation between the positions of differ- 
ent atoms, which is true for gas atoms or dopants in a crystal, so that the 
two-particle probability distribution reduces to a product of single-particle 
distribution functions. To evaluate the integral in fl26l) . we recall that 

= /^^e"*'"^'-^-^), (27) 
fcolr^ — Yji\ J Arc 

where n is a unit vector whose direction is uniformly distributed over sphere. 
Then the integrals of the Gaussian type are readily evaluated: assuming that 
the sample size is much larger than the resonant wavelength, 

koRo > 1, (28) 
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we arrive at 

7coi = 7iN{koRo)-^. (29) 

The factor (/cq-Rq)"^ is the effective sohd angle of the collective forward emis- 
sion of a photon. 

This collective process prevails over the incoherent scattering by individ- 
ual atoms if A^(/co-Ro)~^ ^ 1- Since ~ r/i^g, r] being the atomic density, 
this condition is equivalent to the requirement to have a large number of 
atoms in a cylinder whose length is equal to the sample size in the incident- 
photon direction and the cross-section area of is the order of the wavelength 
squared, i.e. riRok^^ ^ 1. 

The states of the {N — 1, 1} symmetry become populated as well. To the 
first order of the perturbative analysis. 



c, 



{N-1,1} 
I 



it) 



-Si 



Si 



'~1col 



1 - e 



(30) 



To perform the calculations with an accuracy ~ l/y N, it is sufficient to 
make in (IMI) the approximation 



f]^N~'-5,i. 



(31) 



Then 



si = -7iViV 



' EF(r.-r,)--^EE^(r.-r. 



N N 



N 



Ar2 



(=1 i=i 



J dhi J dhjV{ri)Virj)F{ri - r, 
Straightforward but lengthy calculations using (I27|) yield in the limit of 



(32) 



c 



{N-1,1} 
I 



it) 



(l - e"^-'*) . 



(33) 



'N{koRo)^ 

If within the time interval t ^ 7^} we do not detect an emitted photon 
by a perfect detector with 100 % counting efficiency, this implies that the 
atoms have been coherently transferred to the state 



N-l 



(34) 



1=1 



Its normalization coefficient can be found in the limit ^ 1 to be 



^ = I ~ I = / ^ C35^ 

From fl551) we can calculate the probability of such an outcome, namely, that 
no photon is scattered forward during the decay of the fully symmetric state 
and, instead, the new state (!34l) is formed and then decays. The probability 
appears to be small, ~ (/cq-Ro)^^ ^ 1, that is of the order of the body angle 
characteristic for the superradiant forward scattering. Yet for mesoscopic 
samples it may be non- negligible, as argued below. 

Now we may find the explicit coefficients of the expansion \^r) = J2jLi hj\ej) 
for flMl) . To the accuracy ~ 1 / ^/N we find that 




hj = iA koTj - — 51 koFj' . (36) 



The ffist term in the brackets is of the order of koRo, while the second term is 
of the order of koRo/^/N. We can therefore set hj ^ iAkoVj. The radiative 
decay rate of the state is then 

N N 

271 ^ ^ 

E E (korj)(korjv)F(rj - r^v) 



Finally, we obtain 



2A^7i 

knRr 



dhj J dhj, V{r,)V{rj,){koTj){kor,,)F{rj - r,v). (37) 



71 

= MM- ^''^ 

The foregoing analysis has shown that the probability of the phonon 
emision ffist decreases as exp(— 27co«t). For t ^ 7"}, an "afterglow" due 
to photon reabsorption into the states of {A^ — 1, 1} symmetry may occur 
with probability ~ (fco-Ro)~^ 6xp(— 27^^). The entire emission process occurs 
into the forward preferred direction, consistently with the results of |9]. 

Since the cross-section of resonant photon absorption is ~ k^'^, and the 
atomic number density is ~ N/R^, the forward-emission enhancement factor 
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A^/(fco-Ro)^ is of the order of the optical density of the sample. Can this en- 
hancement, N/lkoRoy ^ 1, be consistent with uniform excitation probabil- 
ity over the sample [cf. Eq. f|T71) ]? It can, e.g., if we use an auxiliary strong 
laser field acting on a different transition to cause the the Autler-Townes 
splitting of the optical transition at ujq, thereby reducing the optical density 
of the sample at ujq well below 1. After a low-probability non-coincidence 
event (detection of the signal photon only) that signifies the absorption of the 
probe photon in the sample, the auxiliary field should be rapidly switched 
off, thus restoring the large optical density and allowing the Dicke-state dy- 
namics described above. The need to switch off the Autler-Townes splitting 
at the proper time to observe the enhancement of the single-photon emission 
stresses the relevance of the statement [H] that "timing is everything" . 

Our analysis has provided new, more detailed insights into the buildup 
of the coperativity in space and time in the spontaneous emission process 
triggered by the controlled absorption of a single photon [H] at t = 0. It 
has underscored the dominance of the symmetric Dicke state, but only in the 
long-time and large-sample asymptotic regime. Other cooperative states, of 
lower symmetry, become mixed with the symmetric state as the decay process 
unfolds and add slower emission rate in the forward direction. However, their 
contribution is non-negligible only if the sample is mesoscopic, say koRo ^ 10. 
This analysis thus corroborates the results of Scully et al. [9j in the large- 
sample limit. 
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